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Abstract 

In this note, we study an optimal transportation problem arising in 
density functional theory. We derive an upper bound on the semi-classical 
Hohenberg-Kohn functional derived by Cotar, Friesecke and Kliippelberg 
[6] which can be computed in a straightforward way for a given single 
particle density. This complements a lower bound derived by the afore- 
mentioned authors. We also show that for radially symmetric densities the 
optimal transportation problem arising in the semi-classical Hohenberg- 
Kohn functional can be reduced to a 1-dimensional problem. This yields 
a simple new proof of the explicit solution to the optimal transport prob- 
lem for two particles found in [7J. For more particles, we use our result 
to demonstrate two new qualitative facts: first, that the solution can con- 
centrate on higher dimensional submanifolds and second that the solution 
can be non-unique, even with an additional symmetry constraint imposed. 



1 Introduction 

In this paper, we study a multi-marginal optimal transportation problem arising 
in density functional theory (DFT) in condensed matter physics. Optimal trans- 
portation is the general problem of coupling two (or, in our case, N) probability 
measures together as efficiently as possible, relative to a given cost function c. 
This is a rapidly expanding area of mathematical research, with many diverse 
applications; recent progress is described in the books by Villani [251 126) . 

To precisely formulate our problem, fix a probability measure p on R d (the 
most physically relevant case being d = 3). Let H(p) be the set of all probability 
measures pn on M. dN whose marginals are all p; that is, {TTi)#PN — p for i — 
1, 2, N, where iTi : M. dN — > M. d is the ith canonical projection. We then define 
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Cn(pn) ■= 




and 

E N [p]:= inf C N {p N ) (1) 
p«en(p) 

Readers familiar with optimal transportation will recognize this as the multi- 
marginal Monge-Kantorovich optimal transportation problem, with equal marginals 
p and cost function c{x\, x%, xn) := \ x -x \ ' wmcri we wm refer to as 

the Coulomb cost. 

Density functional theory is a modeling method used by physicists and 
chemists to understand electron correlations. It was originally proposed by 
Hohenberg, Kohn and Sham [H][T5]i see [T7] and [5] for a detailed introduction. 
The Hohenberg-Kohn functional plays a central role here. Imagine a system 
of N electrons, interacting via the Coulomb potential. Given a prescribed sin- 
gle particle density p, this functional returns the minimum energy among all 
n-particle wave functions ip whose single particle density is p. The Hohenberg- 
Kohn functional is given by: 



ft f / N\ f 1 

Fhk[p]-= inf (- / V" \V Xi tp\ 2 dp N (x 1 ,x 2 , ...,x N )+[ n I / , ;dp 2 (x 1 ,x 2 )) ■ 

M>^ P V 2m e J RdN f-^ \2J J R 2d \xi - x 2 \ ' 

Here, h is Planck's constant over 2tt and m e is the mass of an electron. 
The notation ip — >■ p means that p is the single particle density corresponding 
to the wave function tp, and p^ and p 2 are the N and two particle densities 
corresponding to ip, respectively. The first term is the quantum mechanical 
kinetic energy while the second is the Coulomb interaction energy. Numerically, 
this expression is unwieldy for large N, as the complexity of minimizing over the 
space of N particle wave-function grows exponentially in N. It turns out that 
the kinetic energy term can be dealt with relatively easily (see [7] for details); 
one of the major goals of DFT is to approximate the interaction energy of two 
electrons by a function of the single particle density p. 

Two recent papers by Cotar, Friesecke and Kliippelberg [7] [6], as well as 
a paper by Buttazzo, De Pascale and Gori-Giorgi [2], have revealed interest- 
ing connections between this problem and optimal transportation^ Of par- 
ticular interest to us in the present work, the work in [6] showed that in the 
semi-classical limit, h — > 0, the Hohenberg-Kohn functional reduces to E^. In 
particular, contributions from the kinetic energy vanish and the antisymme- 
try property associated with iV-body wave functions reduces to symmetry of 
the measure p^ in the arguments [x\, X2, Xn). In fact, given any measure 
Pn G n(p), we have Cn(pn) = Cn(p~n), where p~N is the symmetrization of pn] 
that is, for all Borel A C R Nd 

1 In fact, these connections were already implicitly present in the physics literature, but 
without rigorous justification |24l| 23 | . 
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ph(A) = — ^ PN{{{x„(l) : X <j{2 )T--X < y{N)) ■ (x 1 ,x 2 ,--,x N ) G A}), 
' creS N 

where Sn is the permutation group on N symbols. Therefore, when computing 
En[p], we can neglect the symmetry condition entirely. 

Optimal transportation problems with two marginals (ie, N = 2) with a 
variety of different cost functions have been studied extensively, including results 
in [7J and [2]on the Coulomb cost; see, for example p] [10] [9] [3], or, for a detailed 
overview and extensive bibliography, see [25] and [25] • However, relatively little 
is known about optimal transportation problems with more than two marginals, 
except in certain special cases pi ] [5][n][3][2l3 ] [ro] [ 2T ] p2] P ^ ] . 

In particular, despite its importance in density functional theory, to this 
point very little is known about the structure of minimizers of Cjv for N > 2, 
or the minimal values E^[p\. The main goal of this note is to contribute to this 
problem. 

This paper features two main contributions. The first is the construction of 
a measure pn € 7r(p); the immediately yields an upper bound on Ejq[p], which 
can be calculated explicitly from the density p. When d — 1, we suspect that 
this measure is optimal and prove this in a special case, although we leave the 
general case open. In higher dimensions, our measure couples every line through 
the origin to itself, using the aforementioned 1-dimensional construction. A 
lower bound on E N [p] was identified in [7J; precisely, E N [p] > {^)E 2 [p\ = 

N< " N 2 — -E^fp]- El For radially symmetric measures when N = 2, Cotar, Friesecke 
and Kluppelbcrg were able to construct the optimal measure explicitly and 
therefore calculate -E^M- Therefore, in this case, we have a lower bound which 
can be calculated and therefore can be compared to the upper bound we derive. 
For one particularly simple p we do this. We note that when N = 2 and p is 
radially symmetric, our construction coincides with the explicit solution in [7j. 
For larger N, it is not generally the optimal measure. 

This may be useful to physicists and chemists, as the explicit construc- 
tion can be used as a starting point for constructing approximations to the 
Hohenberg-Kohn functional. In fact, as is emphasized in [7J, the fundamental 
goal of DFT is to approximate the electron interaction energy (which depends 
on the two particle density) by a functional of the single particle density. Our 
upper bound, which can be explicitly calculated from the single particle density 
p, could be a good starting point for more sophisticated approximations. 

Our second contribution concerns radially symmetric measures. In this case, 
we show that the problem can be reduced to a 1-dimensional optimal transport 
problem with a certain effective cost function derived from the Coulomb cost 

2 In general, this inequality may be strict, as the optimal measure p2 for the two marginal 
problem may not be the two particle density of some TV-particle wave function. As was pointed 
out in [7], this has to do with the representibility problem for two particle density matrices. In 
the semi-classical limit, where wave functions are essentially replaced by probability measures, 
this is the question of the existences of a measure on M. Nd whose 2d-dimensional marginals 
are all equal to p2- 
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(but not equal to it). This has two immediate applications. First, we use this 
result to provide a new proof of the optimality of a construction in 7 and [2] for 
two marginals, which is considerably simpler than the original proof. Secondly, 
we show that this reduction implies that for more marginals the support of 
the optimizer can be more than d-dimensional and be non-unique (even with 
the additional symmetry constraint. This stands in stark contrast to the two 
marginal case, where solutions are known to be unique and concentrated on 
the graphs of functions over x\ and so are essentially d-dimcnsional. Similar 
phenomena, have, however, been observed in multi-marginal problems for other 
cost functions [5] |20j . Let us note that non- uniqueness of minimizers of Cat 
is also ready present in the d — 1 case when N > 2 (see Remark |2. 1.3ft and is 
perhaps not too surprising. On the other hand, I found the fact that there can be 
more than one minimizer which is symmetric in the arguments (xi,X2, ...,xn) 
surprising; from a physical standpoint symmetric measures pn represent the 
semi-classical limits of single particle densities arising from wave functions. 

Aside from these direct applications, this reduction result should be useful 
in deriving numerical methods to evaluate the Hohenberg-Kohn functional, as 
it reduces an optimization problem over measures on M. Nd to a problem over 
measures on M. N 

The idea that the optimizer should come from aligning the radial densi- 
ties optimally, and positioning the particles on fixed spheres to minimize the 
Coulomb interaction is implicit in the physics literature |24j . but has not been 
made rigorous until now. In particular, the explicit construction in the N = 2 
case, though present in [24] and [23], was not proven to be optimal until the 
work in [7] and [2], by a fairly complicated proof, relying on a general theorem 
guaranteeing the existence and uniqueness of Monge type solutions. Our proof 
here, essentially making rigorous the intuition in [23], is much simpler. On the 
other hand, the non-uniqueness we find for larger TV does not seem to have been 
anticipated at all. 

The next section is devoted to the derivation of our upper bound on E^[p\. 
In section 3, we restrict our attention to the radially symmetric case and show 
that the problem can be reduced to an optimal transportation problem with 
1-dimensional marginals. 

I would like to thank Robert McCann for originally pointing out the work 
in [7] to me. I am also grateful to both Robert and Codina Cotar for very 
interesting and useful discussions on this topic, and to Gero Friesecke for a very 
insightful explanation of the form of the semi-classical limit of the Hohcnberg- 
Kohn functional. 

2 An upper bound on En[p\. 

Here we construct a measure pn £ n(p)j yielding an upper bound on E^[p]. We 
do this by coupling every line through the origin to itself, via a 1-dimensional 
coupling described below. 
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2.1 The 1-dimensional problem 

In this subsection, we consider (pj when d = 1. We construct a measure pjv G 
n(p) which we suspect, but do not prove, attains the minimum in (|TJ) . Given a 
probability measure p on M, we divide the real line R into N subintervals, each 
with equal mass. Our measure will then couple together mass from distinct 
intervals in each of the N copies of R. This construction ensures that all of 
the electrons remain isolated from each other and generalizes the 1-dimensional 
construction in [7] in the radially symmetric case when N = 2. This construction 
is essentially the same as the ones in [2] and [23], but it's use in constructing an 
upper bound for non-radial densities is, to the best of my knowledge, new. 
Set r° = — oo. Now define r l recursively by 

r*=inf{f : p(r*-V) < 1}, 

for i — 1, 2, 3, N — 1. Set r N = oo. Assuming p is absolutely continuous with 
respect to Lebesgue measure, this yields N subintervals of R, P = [r* -1 ,r*], 
each with p(P) — jj. Define F % : I 1 — »• P implicitly by 

p(r i -\F i (t)) = p(-oc,t). 

Then, letting p 1 be p restricted to P, we have F^p 1 = p 1 by construction. It 
is also worth noting that each F % is an increasing function, and that F 1 is the 
identity function, = t. 

Now, we can extend F l : R — > R as follows: take the image F l (I^) to be 
(here addition is modulo N, so that, for example, I N+1 = J 1 ). F % : P — > 
i i s defined implicitly by 

p(r i+j - 2 ,F i (t))=p(r j -\t). 

It it then clear that {F l )#p — p. Setting F — (F 1 , F 2 , ....F N ), the measure 
Pn ■= F#P is in II(p). 

Remark 2.1.1. (Local optimality of pn) Note that the measure pn is con- 
centrated on the union of sets of the form P l — P x P +1 x ... x p+ N ~ 1 ! f or 
i = l,2,...iV (here, as above, the indices j in P should be understood modulo 

Ft 2 

N). Note that 9 ^ig x j < for all i ^ j. Now, consider the optimal transporta- 
tion problem on P l , with marginals p 1 , p l+1 , p 1+N ~~ 1 . A theorem of Carlier 
J^j, rediscovered in Jiff) /, demonstrates that the measure pn\p^ = P#Pi above is 
the optimal measure. 

We suspect that p^ is optimal, but do not prove this in general. Below, 
we prove this is the special case where p is uniform measure on [0,1]; a similar 
argument for the N = 3 case was worked out in [2]. 

Theorem 2.1.2. Suppose p is uniform measure on [0, 1]. Then p^ is optimal. 
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Proof. As p is uniform measure on [0, 1], the interval P = [(i — 1)/N, i/N], and 
the map F l (xj) = x 3 + (i + j — l)/N, for Xj £ P . 

Define a continuous function u : [0, 1] — >• M as follows. Set u(0) — 0. For 
x £ P, set 

du -N 2 x ^ A 2 

^ ^ (i - j) 2 + ^ (» - j) 2 

Note that this implies ^ (x) is constant on P , and that ^ (a;) is non- increasing, 
so that u is piecewise linear and concave. 

We will now show that c(x\,X2, ...,Xn) — X^ili u ( x i) attains its minimum 
value at every point in the support spt(piv) of p^. This will imply the desired 
result, as for any measure 7 £ II (p) we have 



. . N N 

j c(x 1 ,X2,...,x N )d~/ = / [c(xi,x 2 ,...,x N ) - S^u{xi) + V]u(a;i)]d7 

N 

> M + V / u(x))dp(x) 



with equality when ■y = pn (here M is the minimum value of c(xi, x 2 , xn) — 

E*=i u Oi))- 

Consider the set 

S = {J{(F^ 1 \x),F°W(x),...,F°W(x))\x€[0,l]}, 

(7 

where the union is over all permutations a on 1, 2, ...,7V. Note that spt(pN) is 
contained in S. It is straightforward to see that c(x\,X2, ijv) — J2iLi u ( x i) 
is constant on S; we now show that any point where this function is minimized 
must belong to S, which will yield the desired result. 

As c = 00 whenever x% — Xj for any i ^ j, the minimum of this function 
must be attained at some point where Xi ^ Xj for all i ^ j, ie, at some point 
where c is smooth. We must therefore have -§^-(xi,x 2 , ...,Xn) — ^{xt) = for 
all i. We will show that this implies (xi, ...,xjy) £ S. 

First of all note that, for (xi,x 2 , --^Xn) £ S, a straightforward calculation 
confirms that -^-{x\,X2, ...,Xjv) — j^( x j) = for all j. 

Now, choose a point (xi, X2, xjv) such that -§§-{x\, X2, xjy ) — ^ (x j ) = 
for all j and assume, without loss of generality, that x\ < x 2 < ... < xjy. By 
the pigeon hole principle, we have that Xi £ P for at least one fixed i. Now, 
note that for fixed tji = Xi, the function 

N 

(yi, -Vi-uVi+i, -,j/jv) i-> c(j/i, ...,j/jv) - 
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is strictly convex on the region 

{{yi,...yi- U yi +1 ,...,y N ) : yj < y j+ iVj = 1, 2, ...N} 

and so there is only one point (yi, ...yi-i, y%+i, y n ) in this region where 
-§§r(x\, X2, Xn) — j^( x j) — 0. We already have this equality when yj = 
(j — i) + Xi for all j ^ i and so we must have Xj = (j — i) + Xi, which means 
(x\,X2, • xn) € S, as desired. 

□ 

Remark 2.1.3. (Non-uniqueness of the optimal measure) It is interesting to 
note that the optimal measure in this case is not unique for N > 2. Indeed, note 
that pjv is not symmetric; it's symmetrization p~N is therefore another optimizer. 
However, p~N is the unique symmetric optimizer. In the next section, we will 
see that in higher dimensions, there can in fact be more than one symmetric 
optimizer. 



2.2 An upper bound in higher dimensions 

We now construct our measure, by essentially coupling each line through the 
origin with the same line in each of the other copies of R d . This coupling is 
done via the 1-dimensional construction above. 

It will be convenient to work in polar-like coordinates, as M. d is the (almost 
disjoint) union of lines through the origin. Note that, neglecting the origin, 
R d = K x P d_1 , where P d_1 denotes n — 1 dimensional projective space. Let 
r e K and 9 E P d_1 represent the angular coordinates. Note that fixing 
corresponds to fixing a line through the origin. Now, expressing the density 
p(r, 9) is these coordinates, we have, by Fubini's theorem, 

l=[ ( f p(r,9)dr)d9 
Jpd-i V J R ) 

and so r ^ p(r, 9) is intcgrablc for almost all 9. For each fixed 9, and m = 
1, 2, N let r i-> F m (r, 9) be the optimal map obtained in the previous section, 
coupling the density r n> p(r, 9) to itself. For m — 1, 2, 3, ...N, define T m (r, 9) = 
(F m (r, 9), 9). 

Proposition 2.2.1. T m pushes p to itself. 

Proof. Note that T m is clearly bijective almost everywhere and it is smooth 
almost everywhere, as F m is smooth on the interior of each interval P. The 
derivative of P m , in (r, 9) coordinates, is (in block form): 



DT m (r,9) = 



D r F m (r,9) D e F m (r,9) 
/ 



The determinant of this upper triangular matrix is D r F m {r, 9) = p ^ Fm ^ r g) g ^ 

(p(T'"'(r > e)) ' as F m { r iG) pushes forward the density p{r,9) to itself, for fixed 9. 
By the change of variables formula, (T m )#(p) = p. 

□ 
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Now set 

p N = (T\T 2 ,T 3 ,...,T N ) # p. (2) 

By the preceding proposition, this yields a measure on (R'')^ whose marginals 
are all p; that is, an element of n(p). We immediately obtain 

Corollary 2.2.2. For the pn defined in equation fjj|) ; we have: 

E N {p) <C N (p N ) = I V, -dp N {x x ,x 2 ,...,x N ) = I V 1 rr dp(x) 

Remark 2.2.3. For measures which are not too spread out in the angular direc- 
tions, we suspect our construction is nearly optimal. Indeed, in the limit where 
the support of p is a single line through the origin, the problem (Qp reduces to 
the 1-dimensional problem addressed in subsection 2.1. On the other hand, our 
construction is far from optimal for measures which are not spread out in the 
radial direction. Consider the limit when p is concentrated at a single radius 
r. In this case, the one dimensional "densities" r i— > p(r,9) are each the sum 
of two Dirac masses; therefore, our construction yields an infinite total energy 
when N > 2 and is thus far from optimal. 

2.3 Example: uniform radial density 

When N = 2 and p is radially symmetric, Cotar, Friesecke and Kliippelberg 
showed that the pn constructed above is in fact the optimizer in (JlJ; that is, 
E 2 [p] = C 2 {p 2 ) CO- They also noted that, for general N, E 2 yields a lower 
bound on E N : E N [p] > ( N 2 )E 2 [p] = N(N ~ 1) E 2 \p}. For radially symmetric 
measure, then, we now have both an upper and lower bound on En which can 
be calculated explicitly. Below, we compare these two bounds for a simple, 
radially symmetric p. 

Proposition 2.3.1. Suppose that p is the radial symmetric density given by 
p(r) = |X[-i,i]. Then E 2 [p] = 1 and C N [p N ] = =f + f + ^(E^i* £), 
where p^ is as in (0). 

Remark 2.3.2. Note that J2n=i < 1 + ln (^ - 1)- Thus , we have 

r , „ , s N N 2 ln(N - 1) 
E N [p] < C N {p) <- + y - 1. 

On the other hand, 

, , /N\ , , N(N - 1) 
E N [p]> (jE i \p] = -^ '-. 

The ratio of the upper to lower bound, then, is roughly logarithmic in N. 
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Proof. For all x, we have one F m (x) in each interval P. Note that, if F rn (x) € P 
and F n (x) 6 P, we have |F m (ir) - F n (x)\ = 211=21. We then have: 

X ^ 1 X , AT 

^ |F m O) - F n (x)| ~~ ^ 2\i-j\ 

It is straightforward to calculate that there are N—l pairs (i, j) where = 1, 
N — 2 pairs where |i — j| = 2, etc, ending in one pair where \i — j\ = N — 1. 
Thus the above is equal to 



TV NN-l N-2 L 



^2|i-7| 2 V 1 2 iV-f 

Mi 

AT, TV AT AT 

= y^- 1 + y- 1 + y- 1 + - + A^T- 1 ) 

= T^ + Y + y + - + iv-T-( iV - 1 » 
= ==!(i + i + i + ... + ^_)-=V-i) 

2 1 2 3 AT-r 2 1 ' 

2 ^ n + 2 2 

n=l 

Integrating the constant function X^m^n |F"»(x)-F"(z) a g ains t uniform measure 
on [—1,1] now yields the desired result. □ 



3 Reduction to a 1-dimensional problem for ra- 
dially symmetric measures 

We show in this section that, for radially symmetric measures, the problem 
can actually be reduced to an optimal transport problem with 1-dimensional 
marginals. Set 

h(n,r 2 ,...,r N ) = inf c(x 1 ,x 2 , —,x N f^ (3) 

1 2:1 |=)-i , \x 2 I =r 2 , . . . , \x N |=rjv 

Note that in this section, we will work in polar coordinates rather than the 
polar-like coordinates used earlier. That is, we will identify R d w [0, 00) x S^ 1 . 
Radial symmetry simply means that the density p(r, 0) — p(r) of the marginal 
p = p(r, 9)drd9 is independent of 9. 

For each fixed r = (n,r 2 , — , rjv) 6 K^, fix 

(xi,x 2 , ...,xjv) = xi(r),x 2 (r), ...,x N (r) 

3 We take c(xi,X2, ...,Xrf) = 1 — " — r nere i although a similar reduction will hold for 
other rotationally invariant costs. 
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attaining the infimum in (J3|); note that for any rotation A G SO(d), 



(A(x x ),A(x 2 ),...A(x N )) 

attains the infimum as well. Take S ri = {x G K 3 : \x\ = r,}. Now let v be Haar 
measure on the rotation group SO(d) and let ^ri,r 2 ,...,r N _ Q^ v ^ where 

G : 50(3) -> 5 n x SV 2 x ... x S rN 

is defined by G{ A) = {A( Xl ), A(x 2 ), A{x N )). Then Y ur2 '-' rN is a probability 
measure on S ri x S r2 x ... x 5 rjv with uniform marginals. 

Then consider the optimal transportation problem with 1-dimensional marginals 
jit := p(r)dr and cost function h; that is, minimize 

/ h(ri,r 2 ,...,r N )d-f (4) 

among measures 7 on whose 1-dimensional marginals are all p. 

Let 7 be an optimizer in ^ (existence is guaranteed by standard results in 
optimal transport theory) and set 

p N (xi,x 2 ,...,x N ) = PN(ri,,r 2 ,...,r N ,9i,62,...,9 N ) 

=: 7(r 1 ,r 2 ,...,r N )®Y ur2 '-' rN (0i,0 2 ,...,d N ). (5) 

Theorem 3.0.1. Assume that p is radially symmetric. Then the pn defined in 
[5| is optimal for (QJ). 

Proof. It is easy to see that Pm has the appropriate marginals. Now, the Kan- 
torovich duality theorem [13] and the optimality of 7 imply the existence of 
functions u%, ...,un : (0, 00) — > M such that h(r±, rpj) — J2i=i u i( r i) > with 
equality when (r%, rpj) G sptfi). We then have, for all X\, x 2 , xn 



N N 

c(xx,...XN) — y]ui(\xi\) > h(\xx\,..., \x N \) -y]uj(\xj\) 

i=l i=l 

> 

with equality when xx,...,xjsr attain the infimum in (|3]), and \xi\,...\xn\ G 
spt(j); that is, when (x±, xn) G spt(pN)- Now, integrating with respect 
to any p~N G n(p), we get 

- N N 

/ c(xi, ...x N )dp~N > / Mi(|sj|)dplv(si)3:2 ) ...Cjv) = y^ l Ui(\xi\)dp(xi) 
Jm Nd ~[ Jm. Nd ~[ Jn d 

Now, noting that we have equality for p]v = Pn and the right hand side is 
independent of plv G n(p) yields the desired result. □ 
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3.1 The two marginal case 



As a consequence of Theorem 13.0.11 we obtain an easy proof of the following 
result, already implicit in |23j (without rigorous justification) and also present 
in [7] and [2] (with more complicated proofs). 

Corollary 3.1.1. When N — 2 and p\ — pi(r) is radially symmetric, the solu- 
tion is unique and is concentrated on the graph of the function x <-¥ — j|j/(|x|), 
where f : (0, oo) i— > (0, oo) is defined implicitly by 

r—OO rf(r) 

/ pi(s)ds = / pi(s)ds (6) 

Jr JO 

Proof. It is easy to see that the minimum in (J3J is uniquely attained at antipo- 
dal points x 2 = and so h(n,r 2 ) = Noting that Q fJ r , > 0, a 
classical result in optimal transportation implies that the unique minimizer for 
the reduced problem ^ is concentrated on the graph of a decreasing function, 
r 2 = f(ri). As there is exactly one decreasing function pushing the measure \i 
to itself (namely ([5])), this implies the desired result. □ 



3.2 The multi- marginal case 

Finally, we use Proposition 13.0.11 to assert two new qualitative facts about the 
N > 3 case. First, we show that the solution may be concentrated on a set 
with dimension greater than d (note that Monge type solutions, and their sym- 
metrizations, would have dimension d). Secondly, the solution may be non- 
unique, even with an additional symmetry condition imposed; that is, there 
can be more than one minimizer in (fTJ) which is symmetric in the arguments 
(xi,x 2 ,:.,x N ). 

Our results in these directions (Propositions 13.2.11 and 13.2. 4[) rely on certain 
assumptions on the structure on the support of the minimizers. We suspect that 
these conditions hold generally, at least when the number of electrons is large. 
Below (Lemma 13.2.2)) we provide an explicit example of a marginal p for which 
optimizer satisfies the conditions in Proposition 13.2.11 Numerical calculations 
in [24] suggest that the conditions in both Propositions 13.2.11 and 13.2.41 hold for 
other marginals (see Remark l3.2.5l b elow) . 

Proposition 3.2.1. Suppose N > 3 and d > 3. Let p be a radially symmetric 
measure, and assume that for some optimizer 7 in ^ the following condition 
holds. 

For all r in some subset I C (0, 00) of positive measure, there is at least one 
point r = (ri, r 2 , ...rjv) £ spt("f), with n = r for which there exist minimizing 
(x\(f), X2if), xisr(r)) in ([3]) such that Xi(r) and Xi(r) are not co-linear, for 
some i. 

Then there exists an optimizer pn in (JlJ whose support is at least (2d — 2)- 
dimensional. 
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Proof. Choose (xi(r), x 2 (r), xjs[{r)) as in the hypothesis; without loss of gen- 
erality, assume x\(r) and x 2 (r) are not co-linear. Now, it is easy to check that 
the set 

{A{x 2 {f)) : A e SO(d), A( Xl (r)) = Xl (f)} = {y G W l : yxi{r) = x 2 {f)-xi{r) and \y\ = \x 2 (r)\} 

is at least (d — 2)-dimensional. Therefore, the set 

{(A(xi(r)),A(x 2 (r)),....,A(x N {r))):AeSO(d), A( Xl {r)) = Xl (r)} 
= {(x 1 (r),A(x 2 (f)),....,A(x N (r))):AESO(d), A(xi(r)) = Xi(f)} 

is (d — 2)-dimensional. 

To summarize, this implies that for each x\ G M. d with radius \x\\ = r G J, 
there is a (d — 2)-dimensional set of points 

{(x 1 ,A{x 2 ),....,A(x N )) : A G SO(d), A( Xl ) = Xl )} 

in the support of the optimizer; the result follows. □ 

We suspect that the condition on the non co-linearity of the optimally cou- 
pled vectors x%, x 2 , Xn is generically true. Below, we prove that this condition 
holds for a specific example marginal p. 

Lemma 3.2.2. Let N > 3 and fix r 1; r 2 , r w € [1 — e, 1 + e]. Suppose 
Xi, x 2 , Xn € argmiri\ x A— ri c(xi, x 2 , x^)- Then, fore sufficiently small, 
and any distinct i,j,k, Xi,Xj,Xk are not co-linear. 

Proof. If Xi,Xj,Xk were co-linear, at least two of them, say Xi and Xj, would 
have to point in the same direction, which case \x{ — Xj\ = \ri — rj\ < 2e. But 
then c(xi,x 2 , ...,xn) > i — - — r > tj-, which is a contradiction for small e as 
other configurations clearly have lower total cost. □ 

Corollary 3.2.3. Assume N > 3 and d > 3. Then there exists a measure p on 
M. d for which there is a measure pn, on l$L Nd , optimal in ([1]), whose support is 
at least (2d — ^-dimensional. 

Proof. Choose an absolutely continuous, radially symmetric marginal p, so that 
the corresponding measure p := p(r)dr is concentrated around r = 1; precisely, 

K{r6[l- £| l + e]})>l-i 

with e as in the previous lemma. Then for any measure 7 G n(/i), and any 
i = 1, 2, ...N, we have 

j{{{ri,r 2 ,...,r N ):n [1 - e, 1 + e]}) = M ([l - e, 1 + ef) < —. 
Therefore, 

7(|J{( r i' r 2,-,r-Ar) :r, ^ [l-e,l + £ ]}) < -. 
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Noting that the complement of Ui=i{( r i' r 2i ••-> r Jv) : r i & [1 — e i 1 + e ]} i s 
[1 - e, 1 + e] N , this implies 7 ([1 - e, 1 + e]^) > \- 

It follows that, for the optimal measure 7 in ©, letting / be the projection of 
[1 — e, l+e] N nspt(^), we have fi(I) > and by absolute continuity, I has positive 
Lebesgue measure. By the preceding lemma, the condition in Proposition 13.2. T1 
is satisfied, and therefore, the result follows. □ 

Proposition 3.2.4. Suppose N > 3 and d > 2. Let the marginal p be radially 
symmetric, and assume that for some optimizer in (j4|) there is at least one point 
f= (n, r 2 , T"3, fiv) G spt{pj) for which the following two conditions hold: 

1. The radii r\, r 2 , Tn a^e distinct. 

2. There exist Xi(r), x 2 (r), x^ir) minimizing ([3]) that span M. d . 

Then there is more than one symmetric minimizer p^ m[7J 

Note that the first condition seems to be generic; it would be surprising 
if the optimizer coupled exclusively points with common radii together. The 
second condition cannot hold if the dimension d is larger than the number of 
electrons N. In fact, it is straightforward, using Lagrange multipliers, to show 
that the minimizing vectors x\{r), X2(r), XAr(r) are linearly dependent, and 
so the first condition cannot hold when d = N either. Nonetheless, I suspect 
the second condition holds generically when d < N. 

Proof. Note that can choose the minimizers (xi(r), x%(f), xj^(r)) in ([3]) to be 
symmetric whenever the radii are distinct; that is, for any permutation a on 
1, 2, ...N, we can assume 

x i (r 1 ,r 2 , ...r N ) = aV(j) (r CT(1 ) , M2)-MiV))- 

Now, the procedure above produces an optimal measure pn, defined by 
([5]); note that we can find another minimizer pjy simply by replacing in the 
procedure leading up to (0 Haar measure v on SO(d) with Haar measure on 
the whole orthogonal group 0{d). Let p and ~p be the symmetrizations of p 
and p, respectively. We need only to verify that p~N ^ pw- To do this, we will 
exhibit a point in the support of p~N which is not in the support of p~N . 

Choose f*= (ri,r2, ...,r]v) € spt(j) with distinct radii and (xi,X2, ...,xn) := 
(xi(r), X2(f), XN{rj) G spt(pN) which span R d . It follows that for any A € 
0(d) such that A £ SO(d), we have 

(A( Xl ), A(x 2 ), ...,A(x N )) ^ (A( Xl ),A(x 2 ),...,A(x N )) 

for all A (^SO{3) (as if (A(xi), A(x 2 ), A(x N j) = {A( Xl ), A(x 2 ), A(x N )), 
then A — A by the span condition. It now follows that 

(A(xi), A(x 2 ), ...,A(x N )) 6 sptij^) C spt{j^) 

Now, it follows from the construction that the only points Z\, z 2 , z?f in the 
support of pn with \zi\ — ri for all i = 1,2,.., N are points of the form 
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(A(xi), A(x2), ■ A(xn)) with A G SO(d). The only points z\, z%, zjv in the 

support of pat with \zi\ — n are points of the form (A(x^ 1 ) (f CT )), A(x CT ( 2 ) (rv))> A(av(jv-)(*v))) 

with ^4 G SO(d), where tr G S^r is a permutation and 

»V = 0V(1))*V(2), •■•' r o-(A')) 

But as the minimizers x 2 (r)..., ^iv(^)) are symmetric, (v4(x CT (i) (r CT )), A(av(2) (fa)), — > A(x a ( N ) (rv))) = 

(A(xi), A(x2)j •••)-A(aJiv))) and so 

(A(xi), A(x 2 ), A(x N )) <£ spt(p N ), 

from which the result follows. 

□ 

Remark 3.2.5. (Numerical evidence in support of the hypotheses) In 

I24L the authors compute what they call co-motion functions. Using spherically 
symmetric densities p on R 3 corresponding to the ground states of the Lithium 
atom (N=3) and Beryllium atom (N=4), they find functions fi, i = 2,...,N, 
such that, in our notation, 7 := (Id, / 2 , //v)|#/i G n(/x). They also numeri- 
cally construct potential functions v(x) and verify numerically that 

N N 

(xi,x 2 , -,x N ) i-> ^2 1 — ~ — 1 - ^2 v ( x *) 

i^j \ xi x l\ i=l 

is minimized at each point on the set {x, f 2 (x), f]y(x)}. This easily implies 
that, in our notation, 7 is optimal in Q. They also compute the minimizing 
angles numerically, and one can readily verify from their calculations that: 

1. For all x, the radii in \x\, |/ 2 (x)|, |/at(x)| are all distinct. 

2. For the Beryllium density with N — 4 7 the vectors x, f%(x), ...fjsr(x) are 
linearly independent. 

Provosition \3.2.7\ then implies that for these densities there are optimal mea- 
sures pn with at least 2d — 2 = 4- dimensional support, while Proposition \3.2.~4\ 
implies that there is more than one symmetric optimal measure for the Beryllium 
density with N = 4. 
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